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O il We propose a scheme to simulate the ID Majorana equation with two Cooper pair boxes coupled to 



a ID superconducting transmission line resonator, where strong coupling limit can be achieved. With 
proper chosen of systematic parameters, we are able to engineer different kinds of interaction, which 
' (— I ' ' is indispensable in simulating the Majorana equation in an enlarged real Hilbert space. Measurement 

of the conserved observable of pseudo-helicity via transmission spectrum of the cavity field can 



^ , verify the simulated Majorana wave function. The measurementresults are experimentally resolvable 

^ i based on our estimation with conservative parameters. 
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o' 



The request for a quantum simulator is because that it allows for the study of quantum systems that 



'sj- \ beyond the reach of classical computers. Meanwhile, it is expected to be more robust against virous im- 

perfections than a quantum computer ||1|]. One of the successful cases is the quantum simulation of the 



Dirac equation, which combines quantum mechanics and special relativity 



m. 



In certain regime, electrons 



KA ' in graphene may behavior as Dirac fermions, which has recently raised strong interest in condensed-matter 

C^ . physics ISO. Meanwhile, it is proposed that ultracold atoms in an optical lattice can be used to simulate 

relativistic Dirac fermions [4]. Compared with graphene, the atomic simulation may offer more controlla- 
bility. Alternatively, simulation of the Dirac equation is also proposed in cold atoms ||5|] with light-induced 
gauge potential laLZl]. Meanwhile, quantum simulation of the Dirac equation in trapped ions system is also 



proposed ISf] and been experimentally confirmed 



m. 



One of the greatest success of the Dirac equation is that it predicts the existence of antiparticle of elec- 
tron. In viewing of the success of the Dirac Equation, Majorana inquired whether it is possible that a particle 
to be its own antiparticle. As a result, he found an equation that such particle would satisfy, i.e., the Majo- 
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m. 



rana equation (ME) lHOll . Recently, it is proposed that ME can be simulated with trapped ions llllll . which 
is developed in an enlarged Hilbert space: an ID Majorana equation is transformed to a 3D Dirac equation 
with dimensional reduction, i.e., momenta in y and z directions are zero. The simulation of the ME is not 
straightforwardly as it is non-Hermitian, i.e., one needs to implement the complex conjugation of Majorana 
wave function. Therefore, it needs new toolbox to access the unphysical operations in the simulation. 

Superconducting system is generally regarded as one of the most promising candidates for physical im- 



13 1 ■ Recently, quan- 



plementation of qubit that can support scalable quantum information processing 11121 
tum simulators using superconducting quantum circuits have attracted much attention II14I1 . With mature 
microchip fabrication techniques, it can be used to simulate quantum system in a flexible way. Here, we 
consider the simulation of ME with two Cooper pair boxes (CPB) capacitively coupled to an ID transmis- 



sion hne resonator (TLR) 111511 . where elements are mesoscopic circuit excitations. Therefore, it provides 
an interesting example of simulating physics associated with relativistic ME in a mesoscopic circuit. As 
the simulation in Ref. [11] is based on complex-to-real map, which transforms a Majorana equation into 
a higher dimensional Dirac equation, for 2D or 3D ME simulation, one needs more qubit (more than 3) to 
simulate the transformed Hamiltonian. Unfortunately, in trapped-ion systems, the individual addressing is 
difficult for large arrays (A^ > 3) jlol . Comparing with the trapped ion simulation, a distinct feature of the 
present proposal is that the combination of individual addressing with a many CPB setup is feasible. 
A ME reads 

ih'y'^d^'ip = mcipc, (1) 

where c is the speed of light, m is mass, 7^ with ;U=1, 2, 3, and 4 are the Dirac matrices, and tpc = i^'^ICi/j is 
a charge conjugation field with /C being a complex conjugation operator. For ID case, tpc = Cazip*, where 
C is a unitary matrix satisfying Cj'^ = — 7C with 7 = iaxCTz- Then we can choose C = ioy in a suitable 
basis, i.e., ipc = icTyUzJCip = iay(Tzip*. Then, ME in ID reads 

ihdti^ = {caxPx — iiTT-c ay)C)il^. (2) 

where ^ is a two-component complex spinor, px = —ihdx is the momentum operator in x direction. 

To make the nonphysical complex conjugation ip — )• JCtp = ij)* implementable, one may map the two- 



component complex spinor into a four-component real spinor lllin . that is, 

/ Vi \ 1 / V + v-* \ 

V' = G C2 ^ * = - G M4. (3) 

After this mapping, one can unify all the antiunitary or unitary operators and complex field V' in an enlarged 



Hilbert space. Then, in the enlarged real Hilbert space, Eq. Q reads 

ihdt-if = [c(E (g) a^)px - mc^a^ (g) ay]"^, (4) 

where ^ = (-^i "02 V'l "02 ) ^ ^4 ^^'^ ^^ ^'^^ become an implementable Hamiltonian equation. In 
Eq. ^, the four components of the spinor is nonlinearly coupled and can not be separated, so it is a 1 + 1 
Dimension Dirac equation with an irreducible 4 Dimensional Hilbert space notation. In the following, we 
will simulate the dynamics of ME in 1+1 Dimension by engineering this Hamiltonian. Meanwhile, it is 
well known that Majorana equation in 3D conserves the observable of helicity. In ID, helicity reduces 
to so-called pseudo-helicity S = axPx, and it is still conserved in ID ME but not in ID Dirac equation 
ihdt'^ D = Hd'^d with H£) = caxPx + mc^az as [S,///)] ^ 0. Therefore, measurement of the pseudo- 
helicity can illustrate the Majorana wave function, and thus to verify our dynamical simulation. Mapping 
into the real Hilbert space, the pseudo-helicity S reads 

t = M^axPxM = (E(gax-cry0 ax)Px (5) 

where E is identity matrix. 

We now turn to our circuit QED simulation. The qubit considered here is the superconducting Coopper 
pair box (CPB) consisting of superconducting island where two Josephson junctions with capacitance Cj 
and Josephson energy Ej are configured into a loop geometry, which is pierced by an applied external 
magnetic flux $. On condition that Josephson energy is much small than charging energy Ec = e^/2Cs 
(Cs = Cg + 2Cj), restricting the induced charge Ng = CgV^'^/{2e) to the range of Ng G [0, 1], only a pair 
of adjacent charge state on the inland are relevant. Then, the CPB is a simple two-level system described 



byllJ 



Eel - Ej _ 

where ax, (Jz are thepauli matrices in the basis of {|0), |1)}, E^i = 4:Ec{l — 2Ng) is the electrostatic energy 
and Ej = 2Ej cos(7r<l>/<I>o) is the effective Josephson energy with <l>o being the flux quanta. From Eq. (O, 
one can see that it is possible and convenient to control the qubit by the applied gate voltage V^'^ and the 



pierced flux $ 111311 . Therefore, the qubit splitting energy can be tunable by the external magnetic flux even 
with fix gate voltage, e.g., at the degeneracy point Ng = 1/2. 

In circuit QED, the CPB is capacitively coupled to the center conductor via a capacitance Cg, at the cav- 
ity mode's antinode with a maximum voltage. Then, besides the dc control voltage V^'^, the gate voltage on 
the CPB also include an ac part from the oscillation of the cavity mode. Taking both parts into consideration 
lead Eq. Q to 

Hi„t = -2E, (l - 2A^f ) az - ^a, + hwra^a -hg(a + a^) (l - 2iVf - a,) , (7) 



where a^ and a is the creation and annihilation operator of the cavity mode and Wr is its frequency, g is the 



coupUng strength with g/2TT £ [5.8, 100] MHz ilTP . Then, denoting that {| J,), | t)} are the ground state and 
the excited state of the first two terms of Hamiltonian dV), respectively. In this new basis, at the degeneracy 
poin t, within the rotating-wave approximation, HamiUonian ^ reduces to the Jaynes-Cummings model 



I117D 



Hjc = hwrO^a + —CTz + hg{a^a- + aa+) (8) 

where uj = Ej /h and az is pauli matrices in the new basis. 

In a circuit QED system, a microwave drive of frequency w^ can be capacitively coupled to the resonator, 
which can be in the form of 



Hd{t) = Ht) 



a^f,-i{^dt-^) + ae^("'''*- 



(9) 



where e{t) and (p is the amplitude and the initial phase of the microwave pulse, respectively. All of the local 
operation on the qubit are rely on Wd, £{t) and phase (/), which have been experimentally achieved [15]. For 
the large drive amplitude compare to the vacuum fluctuation of the resonator, the microwave pulse can be 
treated as a classical field. Make a unitary transformation U = exp(aa^ — a* a) on the total Hamiltonian 
consists of (IHD and IM leads to lllSn 



Hdjci = hwrO^a + — (Tz + hg 



a) + a*)a- + H.c. 



(10) 



On condition that the drive amplitude are independent of time, — qA = eeyi'p[—i{wdt — (f))] with A = 
ojr — iOd- In the rotating frame at frequency Wd, we can rewrite (fTO|) as 

Hdjc2 = hAa^a + hg{a)a. + aa+) + ^-a, + ^{a_e-'^ + a+e'^), (11) 

where Qd = 2(7e/A and 6' = uj — ujd- 

Working in the eigenbasis of the last two terms of the above Hamiltonian, the Hamiltonian pT]) becomes 

Hdjc3 = hAa^a + —-a^ + ^[ae"*'^(cos(6l')cT^ - sm{e')a^ + cr+ - o-_) + H.c], (12) 



ISD reads 118, 



m 



where 6' = arctan((5'/ilrf) and Q = \/^d + ^ ■ ^^ ^^^ interaction picture, the interaction Hamiltonian of 



Hj = !^[ae-'^'e-'^{az + cr+e*™ - a_e-*™) + H.c], (13) 

which is key to achieve our effective Hamiltonian for simulating ME and we have 6' = 2n7:Q,d with n being 
an integer. From this microwave driven Hamiltonian, one can engineer different type of interactions. For 
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FIG. 1: Level structure of our scheme, wi, u!2 are level separation of the two CPBs, the frequency of the four 
microwave pulses are ujdi , ^d2, ^dz and ujd4„ and w^i and ujr2 are the frequency of the cavity modes. 



A = ri, one obtains 



H', = ^ae-''f>{a+ + a.e"^^* - a^e-^*™) + H.c. 



(14) 



In rotating-wave approximation, one can neglect the fast oscillating terms, then the above Hamiltonian 
reduces to 



hg 



Hi = -^{aa+e **" + oV_e 



i(l> 



,t. 



P«<A^ 



(15) 



Similarly, for A = — ^2, one obtains 



Ho 



fig. 



{aa.e-'^ + a}a+e'^) 



i<f>\ 



(16) 



Meanwhile, for strong driven Q, ^ A one obtains 

hg 



//, = !|(ae-(At+<^)+ate^(^*+*))^, 



(17) 



To simulate ID ME in circuit QED, we use two CPBs with level separation uji and ti;2, four microwave 
pulses of frequency L0di,i0d2, Wds and com with initial phase (pi, 02, (ps and 04, and two cavity modes of 
frequency cj^i and u;r2- The level structure and the frequency of the driven fields are shown in Fig. 1. As 
explicitly shown in Eq. ^, there are two components: the first one is the kinetic term of CPB 2 cp^i^^Ci^ax-) 
and the second one is the exchange coupling of the two CPBs. As the kinetic term is the combination of 
Hamiltonian (03] ) and (IT6l ). it can be generated by two detuned driving microwave pulses of frequencies 
ujdi and u}(i2 with detuning Ai = ujr2 — ^di = ^i and A2 = ujr2 — ^d2 = —^2, and the initial phases 




FIG. 2: (color online) Expectation values of pseudohelicity S = a^Px in Dirac equation(red line) and Majorana 
equation(blue Hne). 



01 = 7r/2 and (j)2 = —tt/2. In typical experiments, we can choose a;r2 = 10 GHz, w^ii = 9.9 GHz and 
Wd2 = 10.1 GHz, thus it is readily to get J^i = il2 = 100 MHz. Then, we get the combined kinetic term 
Hamiltonian ihg{a^ — a)(E (Jx)/'2, for CPB 2 with level separation of ui2 = 9.7 GHz. Meanwhile, we 
on purposely chosen two CPBs with much different splitting energy with wi = 4.4 GHz, and thus ojdi and 
0Jd2 can not generate similar kinetic term on CPB 1 on the condition that 5i = uJdi — ^ri = 4.9 GHz 
with ujri = 5 GHz, which is much larger than {Ai, IA2I} = 100 MHz. The second term is generated 
by virtual excitation of the TLR, where we need same detuning between the driving microwave field and 
the two CPB. As we choose the different CPB splitting energy, we also need two additional strong driving 
microwave field: {^73,^24} » A = ujri — ujds = ujr2 — ^dA = 500 MHz. Using effective Hamiltonian 
(ITT] ). with initial phase (/ig = ^4 = 0, we obtain, -^ (cr^ (g" E + E (g) az){a^e'''^^ + ae"*'^*). When we choose 
the microwave pulse ujds = 4.5 GHz, and Ud^ = 9.5 GHz, respectively. The cross talk of the two coupling 
channels on CPB2, kinetic and exchange coupling, can be eliminated via rotating-wave approximation 
because of the chosen parameters A S> Ai. Meanwhile, the driven of CPB 1 by Udi is neglecting small 
while 62 = 0Jd4: — ^ri S> A. Then, in the interaction picture, the complete Hamiltonian is 



H = ^{a^(g)E + E(g) a^){a^e''^- + ae" 



+ yi(a^-a)(E(go-x). 



(18) 



Rotating on the CPBl and CPB2, Hamiltonian ([T8j) can be rotated to 
H 



hg_ 
2 



(a^ E + E <Jy){a)e'^^ + ae"^^*) + ^i{a) - a)(E «> a^), 



(19) 



the effective Hamiltonian of which recover the ID ME Hamiltonian in Eq. (HJl with the identification of 



'm'f}hJr2 , t ^ 

n/ — (a' -a), 



h 



2m'Wr2 ' 



mc 



2A' 



(20) 



where m' is the mass of inductance. 



As mentioned above, to verify our simulation, we need to measure the Pseudo-helicity, which is con- 
served in ME. In Fig. 2, we numerically calculated the expectation value of pseudo-hlicity S = axPx in 
both Dirac equation and ME. Obviously, the expectation value oscillate versus time in 1 + 1 Dirac equation 
(the red hne) but not in 1 + 1 ME (the blue line). In our simulation, the initially spinor state is chosen as 

i'(k,t=o) = (i^e-(°'=)' I M , where a = ^^J^ = 0.16 x IQ-i^ with m' = 3.3127 x lO'^^ ^^^ 

ujr2 = 10 GHz, and the FWHM in k space is A A; = ^. For convenience of comparison, we set the the 
third and fourth component of the spinor state initially to be zero in ME. 

We now move to measure the Pseudo-helicity in circuit QED. We first deal with the kinetic part, which 
can be measured as following: 1) a state-dependent operator U2 = exp(— iA;(E ® ay) ® Pxf^) acting on 
CPB 2, which can be generated by two microwave pulses with initial phases being vr; 2) a a^ measurement 
of CPB 2, which can be measured by microwave irradiation of the cavity and then probing the transition 
frequency if the qubit state is in excited state or not B 



ofini] 



1711 . The above two steps equal to the measurement 



F{k) = ul{E (g) a^)U2 = cos{kpx){E (g> a^) - sin(A;p^)(E a^), (21) 

where k is determined by the probing time 191]. As ^(^(^)) L_n = —{(^ ® (^x) ® Px), the kinetic term of 
Eq. (10) can be measured by probing the initial slope of the observable F{k) |9]. Similarly, the second spin 
exchange interaction term in Eq. (10), can be measured by a unitary transformation Ui = ex.p{—ik{ax 
IE) (8) Pa/2) on CPB 1, and measure the spin correlation az CS) ax, for the former operation, an additional tt/2 

= 2{{ay (g) ax) <^Px)- 

Combining both the kinetic and spin exchange terms, we achieve our final goal of measuring the Pseudo- 
helicity in Eq. (5). 

Although the absolute value of Pseudo-helicity is small, it should not hamper our measurement since 
what we directly need to measure is a^. In circuit QED, high-fidelity quantum non-demolition measure- 



pulse is needed. Then probing the initial slope as before, we obtain ^''^j^^ 



ments of az is now an experimental routine Ill7ll . In the follwing, we highlight the measurement process in 
our scheme with input-output process. In our model, if the cavity contains two qubits, the intracavity field 
come from bin will acquire some nontrivial dynamics which then compel the external field aout to have a 
time dependence different from the free field d yna mics (cf. Fig. 3). Then there is a qubit-state-dependent 



phase shift between the input and output fields 112011 . 

As we all know that there is an associated loss mechanism in all physical process, particularly in mea- 
surement scheme, the loss strongly influence the precise quantum non-demolition measurement, such as 
regulating and controlling, so the dispersive regime should be considered in measurement mechanics. In- 
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FIG. 3: A schematic representation of the cavity field, input field and output field for a two-side leady cavity 



sight into the dispersive regime between the CPB 1 and the intracavity can be obtained by the unitary 
transformation 



Ui = exp 



— (ac7+-aV_ 



(22) 



where large detuning A' = Uri — uJi ^ g. Applying transformation (|22] | to the Hamiltonian ([8j and ^ to 
second order in g (neglecting the damping for the moment), for CPB 1, we can get 

Hsysi = Ui{Hjc + Hd)U\ 

, , ^ h hg^ I hg^ hge ^ , \ \ ,^^^ 

ss ^(wri -a;rf3)aia+ -(wi -WdsJCTz + -— a'aa^ + -— CT2 + -— aa; + fe(ai +a^ (23) 



A' 



2A' 



in a rotating frame at the drive frequency w^s. 

We can consider the resonator as a two-side leaky cavity with equal rates, as shown in Fig. 3. The 
relationship between the input and the output mode can be written as aout{t) = yjlia{t) — ain{t), where 
aout and am is the output mode and input mode at the output port res pec tively, a is the intracavity mode, 



and cavity lifetime 1/ k is about 160ns in a typical circuit QED system jlSl ]. The quantum voltage is related 
to the current carried by the TLR by /(x, t) = ^JJV{x, t). For the moment, it is more convenient to have 
stationary rather than traveling quantum voltage. For this mode, we have following equation of motion for 
the CPB 1 

da{t) 



dt 



--[a{t),Hsysi] - Ka{t) + \^ain{t) + ^/nbinit) 

-i[{iOri - ^d3)a + e + xiacTz] - na{t) + ^/Kain{t) + \/K6j„(t), 



(24) 



where xi = 5^ /A' and bin is the input mode at the input port of the resonator. Applying Fourier transfor- 
mation on the above expression and rearranging, we can obtain a spectra function between the input and the 
output mode 



K[ain{uj)+bin{i0)]-l2TT^£5{uj) 

aout{^) = -—, \ ^ ain{w) 



(25) 



9 

where 5{uj) is the Dirac function. Obviously, the Uz in the denominator is a formal indication that the output 
spectra depend on the qubit state {az = ±1). In the case u = ujri — w^a, we can get JllH 

aoutiWrl - Wds) = -— =F -— ain[U}rl - Wrfs). (26) 

K±lXl K±Xl 

We can assume that ain{t) is always the vacuum because of the very small reflection and backscatter of 
the resonator. Furthermore, the two input modes are assumed to be independent and we have neglected the 
two-photo process. Therefore, ain{w) has no contribution to the normally ordered moment aout{w), and 
thus 

\^Zut\^r\ - Wd3))AAl = ( ■—. ) = (e 1 bjnl^rl " ^d?,))M\. (27) 

\ «^ ± «Xl I J^\ 



where (•);v' signify a normally ordered moment, Q^ is the phase shift dependent on the qubit state between 
the input field bin and the output field Oom* [tan(6'i^) = T-^] and 6j„ is a scaled quantity of &j„. When the 
cavity resonance frequency and drive frequency is confirmed, the phase shift dependent on the CPBl qubit 
state is about i^^Mf and the corresponding frequency interval is 29.6 MHz for our scheme, which is readily 



resolvable experimentally lIlSll . 



Similarly, analogical transformation can be applied on CPB 2, we can obtain 

k^l>ut{^'))M2 = /^%^\ = {e^HU-^'))u2, (28) 

\K±IX2 / _sf2 

where uj' = w^i + uJd2 + ^^4 — ujr2, tan (6*^) = -^^ and X2 = 5^ /A" with A" = u;r2 — uj2^ g. For the 
CPB2, 02 = ±^ and the frequency interval is 59.2 MHz. 

In summary, we propose to simulate ME with two CPB capacitively coupled to the TLR in circuit 
QED system where we are able to engineer different kinds of interaction which constructing the wanted 
Hamiltonian for ME. The conserved observable pseudo-helicity is measured via the input-output process, 
estimation based on conservative parameters shows that it is experimentally resolvable. 
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